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ABSTRACT 

This paper reports an experimental investigation of 
bandwidth enhancement of dielectric resonator antennas 
(DRA) using parasitic elements. Substantial bandwidth 
enhancement for the HE 11<S mode of the stacked geometry 
and for the HE 13<5 mode of the coplanar collinear geome- 
try have been demonstrated. Excellent radiation pat- 
terns for the HE 115 mode were also recorded. 

INTRODUCTION 

A dielectric resonator placed over a ground plane can 
serve as an effective radiator [1-2]. Compared to 
printed antennas, the DRA has lower ohmic losses and 
higher radiation efficiency particularly at high fre- 
quencies. However, because of high Q factor, the DRA 
has very narrow bandwidth which severely limits its 
usefulness as an antenna. It has been demonstrated that 
significant bandwidth enhancement can be achieved for 
microstrip antennas by placing closely spaced parasitic 
elements on both sides of or directly above a driven 
patch [3]. It appears that parasitic elements can also 
be used to enhance the bandwidth of a DRA. This paper 
is concerned with an experimental investigation of the 
effects of parasitic elements on dielectric resonator 
antennas . 


THE EXPERIMENT 

The single DRA as shown in Fig. 1(a) is placed on a 
ground plane and aperture coupled to a notched coplanar 
waveguide (CPW) through a slot. The slot is positioned 
along the y-axis (0=90°) so the HE m and the HE 13f modes 
can be excited [4]. Bond wires are used to keep the CPW 
ground planes at equal potential and suppress the 
slotline even mode. The DRA of diameter D and height H 
of 6.223 mm and 2.489 mm respectively is constructed 
from ZrSnTiO,, material of relative permittivity e r2 — 
36.0, and the CPW feed is fabricated from 20 mil (0.508 
mm) RT-5880 Duroid substrate of e rl =2.2. In the ex- 
periment, identical dielectric resonators were placed 
directly above and on both sides of the driven DRA. 


© 1993 IEEE. Reprinted, with permission, from Proceedings of IEEE Antennas and Propagation 
Society International Symposium, Ann Arbor, Michigan, June 1993, pp. 1500-1503. 
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RESULTS AND DISCUSSIONS 

The measured return losses for the single DRA is shown 
in Fig. 2(a). Results indicate that with proper aper- 
ture lengths, either the HE 1U or the HE 134 mode can be 
strongly excited. The aperture length was found to have 
a strong impact on the bandwidth of the DRA. The HE 114 
mode excited with an aperture Lj = 0.40 cm has a 2:1 
VSWR bandwidth of 8.8 %, while the HE I34 excited with an 
aperture L* = 0.68 cm has a bandwidth of only 2.5 %.The 
return losses are less than -10 dB for both modes indi- 
cating strong coupling. Figure 3(a) shows the measured 
patterns for the HE 114 which has a broader main lobe in 
the 0=90° plane than in the 0=0° plane. The pattern is 
asymmetrical along the 0=0° plane due to interference 
from the test fixture. 

Figure 1(b) shows the collinear geomtry where two 
identical dielectric resonators are proximity coupled 
to a driven DRA along the 0=0° plane with one on each 
side of the DRA. The separation between resonators is 1 
mm. It was observed that placing the parasitic dielec- 
tric along the 0=0° plane produces strong coupling with 
the HE 134 mode, while along the 0=90° plane produces 
strong coupling with the HE n4 mode. The strong coupling 
is caused by high field concentrations of the HE 134 and 
HE 114 in the 0=0° and 0=90° plane respectively. Figure 
2(b) shows the measured return losses which are less 
than -18 dB for the HE 134 mode with 0=o° orientation and 
about -50 dB for the HE 114 with 0=90° orientation. The 
parasitic elements has increased the 2:1 VSWR bandwidth 
of the HE 134 from 2.5 % to 3.7 % . The measured patterns 
for the HE n<5 mode as shown in Fig. 3(b) exhibits 
excellent broadside characteristics. 

Figure 1(c) shows the stacked geomtry. The overlaying 
parasitic resonator was found to couple strongly to the 
HE U4 mode where the electric field is maximum near the 
top surface of the driven DRA [6]. As indicated in Fig. 
2(c), the measured return loss is less than -20 dB for 
the HE 114 mode with the higher order mode sup- pressed. 
The parasitic resonator increases the 2:1 VSWR band- 
width of the HE 114 mode to over 5.3 %. As with the 
previous case, excellent radiation patterns have been 
obtained for the HE 114 mode. These patterns are dis- 
played in Fig. 3(c). 


CONCLUSIONS 

Substantial enhancement in bandwidth has been demon- 
strated for the HEn and the HE 13 modes with the para- 
sitic resonators placed directly above and on both 
sides of the DRA respectively. Because of excellent 
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return loss and radiation characteristics, the DRA 
should be excited in HE n mode when used as radiating 
elements in an array. 
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Fig. 1 Schematic Illustrating a Grounded Dielectric 
Resonator Aperture Coupled to a Notched CPW 
feed: (a) Single , (b) multiple and (c) 

Stacked Dielectric Resonator Antenna. 
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( a ) 



(b) 



Fig. 2 Measured Return Losses for: (a) Single, (b) 
Multiple and (c) Stacked Dielectric Resona- 
tor Antenna. 



(a) 



Fig. 3 Measured Radiation Patterns for the HE U * Mode 
for: (a) Single, (b) Multiple and (c) Stacked 

Dielectric resonator Antenna (f= 8.2 GHz). 
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Abstract 

A space power amplifier composed of active linearly 
tapered slot antennas (LTSAs) has been demonstrated and 
shown to have a gain of 30 dB at 20 GHz. In each of the 
antenna elements, a GaAs monolithic microwave integrated 
cir-cuit (MMIC) three-stage power amplifier is integrated with 
two LTSAs. The LTSA and the MMIC power amplifier have 
a gain of 11 dB and power added efficiency of 1 4 percent 
respectively. The design is suitable for constructing a large 
array using monolithic integration techniques. 

Introduction 

The power output as well as the dynamic range of 
microwave solid state devices decreases as the frequency of 
operation increases. Hence to obtain high power at fundamen- 
tal frequencies of several tens of GHz, the output from all the 
devices have to be combined using either conventional power 
combiners or quasi-optical power combiners. 1 In the case of 
conventional power combiners, the combining is done with 
Wilkinson, radial line, and hybrid coupled networks. In the 
case of quasi-optical combiners, oscillators constructed with 
IMPATT diodes, 2 Gunn diodes, 3 MESFETs, 4 or HEMTs 5 are 
integrated with microstrip patch antennas or linearly tapered 
slot antennas (LTSAs) 6 to form an active antenna array which 
combines power radiatively in free space. The advantages of 
quasi-optical power combiners over conventional power 
combiners are higher combining efficiency because of lower 
conductor losses and larger dimensional tolerances with the 
absence of resonance modes. In addition both antennas and 
devices can be integrated on a single semiconductor wafer thus 
simplifying the array construction. The disadvantage of 
oscillator based quasi-optical combiners is that the individual 
oscillators have to be phase locked to a reference source or the 
active array has to be placed in a Fabry-Perot resonator to 
produce coherent radiation. 

Another way to obtain high power is to construct a 
spatial amplifier. One such scheme is the grid amplifier. 7 Each 
unit cell of the grid amplifier consists of a pair of packaged 
GaAs MESFETs with the sources connected together to form a 
differential amplifier and with the gate and drain terminals 
extending radially to form a pair of orthogonal strip antennas. 
Radiation to and from a planar array of identical unit cells is 
quasi-optically coupled by a vertically or a horizontally polar- 
ized beam respectively. 


This paper presents for the first time a spatial ampli- 
fier with GaAs monolithic microwave integrated circuit 
(MMIC) multi-stage power amplifiers. In this approach an 
array of active antenna modules constructed from nonplanar 
LTSAs 8 and GaAs MMIC amplifiers receives signals at lower 
power, and after amplification re-radiates signals into free 
space. The two advantages of the spatial amplifier over the 
spatial oscillator are that only a single stable lower power 
source is required (thus greatly simplifying the combiner 
construction) and that the amplifiers can be individually 
optimized. Figure 1 schematically illustrates a possible 
arrangement for space power amplification. 

Active Antenna Module Characteristics 

The experimental three-element array module is 
showns chematically in Figure 2. The array elements are 
constructed by integrating a GaAs MMIC multi-stage power 
amplifier between two nonplanar LTSAs. 

Linearly Ta pered Slot Antenna 

The feed system of the nonplanar LTSA consists of a 
conventional microstrip with the ground plane tapered to form 
a balanced microstrip. The strip conductors of the balanced 
microstrip are gradually flared with respect to the antenna axis 
to form the nonplanar LTSA. The design of the non-planar 
LTSA and its characteristics are reported in Ref. 8. The 
antenna is fabricated on a 0.02 inch thick RT-5880 Duroid 
substrate. The measured gain of the LTSA is about 1 1 dB at the 
center frequency of 20 GHz. The LTSA has a return loss Sj j 
of 10 dB (2:1 VSWR) over a bandwidth extending from 10 to 
30 GHz. 

GaAs MMIC Multi-Stage Power Amplifier 

The GaAs MMIC three-stage power amplifier was 
designed and fabricated by Texas Instruments for NASA Lewis 
Research Center. 9 A photograph of the amplifier chip is shown 
in Figure 3. The amplifier is constructed on a GaAs substrate 
with an active layer doping level of 2.5x10* 7 cm' . The gate 
widths of the three stages are 1 .2, 2.4, and 6.0 mm, respec- 
tively, and the gate length is 0.5 pm in all the stages. The chip 
size is about 4.0 by 3.0 by 0.1 mm. The bias network is 
incorporated on the chip. The drain voltage V d and 
current I d is 6.3 V and 1.9 A, respectively. The gate voltage 
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Vg is -0.6 V. The measured gain (S 2 j) on a HP 851 OB ANA 
with a 40 dB coaxial attenuator on the drain side of a typical 
amplifier is shown in Figure 4. The gain is greater than 10 dB 
over the frequency range of 18 to 21 GHz. The saturated output 
power measured on a Pacific Instruments scalar network ana- 
lyzer at 20 GHz is about 1 .8 W with a gain of 10 dB and power 
added efficiency of 14 percent. 9 

Experimental Results and Discussions 

A simple measurement procedure has been developed 
to estimate the gain of the space amplifier. This procedure in- 
volves the LTSAs at the input terminals are space fed from a 
single horn antenna while those at the output terminals radi-ate 
into free space. The free space radiation is picked up by a 
second horn antenna which is placed at a far field distance 
from the array. The ratio of the measured received power with 
and without bias to the MMIC amplifiers provides an estimate 
of the gain of the space amplifier. In the setup, the two horn 
antennas are orthogonally polarized but the LTSAs are oriented 
to have the same polarization as their respective horn antennas, 
thus good isolation between the transmitting and the receiving 
horn antennas is established. Also, for comparison purposes, a 
single LTSA was tested as a receive antenna. The H-plane 
pattern is shown in Figure 5 which exhibits a power gain of 6.7 
dB with the MMIC amplifier turned ON. 

Three-Element Array Module 

In this experiment, the LTSAs at the amplifier input 
and output are oriented with the H and E vectors of the re- 
ceiving hom, respectively as shown in Figure 6(a). This 
arrange-ment allows the hom to excite the three LTSAs with 
equal amplitude. The measured radiation pattern is shown in 
Figure 7(a) with the amplifiers turned ON and OFF, respec- 
tively. The gain increases by as much as 30 dB when the 
amplifiers are turned ON which is in good agreement withthe 
measured gain of the amplifiers. 

A second experiment, as shown in Figure 6(b), is 
carried out with the LTSAs at the input and output oriented 
with the H and E vectors of the receiving hom, respectively. 
The measured radiation pattern is shown in Figure 7(b). In this 
arrangement the gain increases by 25 dB when the amplifiers 
are turned ON. The gain is lower in this case because the 
LTSAs on either side of the center element are excited with a 
lower amplitude due to the amplitude taper of the electric field 
distribution of the transmitting hom. The experimental three 
element LTSA MMIC array module is shown in Figure 8. 

Conclusions 

A space power amplifier composed of active LTSA 
antennas has been demonstrated and shown to have a gain of 


30 dB at 20 GHz. In each of the antenna elements, a MMIC 
three-stage power amplifier is integrated with two LTSAs. 

The GaAs MMIC power amplifier and the LTSA have a power 
added efficiency of 14 percent and a gain of 1 1 dB, respec- 
tively. The design is suitable for constructing a large array 
using monolithic integration techniques. 
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Figure 1 . — Schematic illustrating a possible arrangement for 
space amplification. 
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Figure 2. — Schematic illustrating the three-element array 
module. (X 0 : free space wavelength.) 




Figure 4.— Typical measured gain of MMIC amplifier. 



Figure 5, — The measured H-plane radiation pattern of a single 
LTSA with and without the amplifier. 
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(a) H-plane. 



(b) E- plane. 

Figure 6. — LTSA orientation in the three-element array module 
for gain measurement. 




(b) E- plane. 
Angle, deg 


Figure 7.— The measured radiation pattern of the hom antenna 
showing space power amplification. 



Figure 8. — The experimental three element LTSA MM 1C array 
module. 
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ABSTRACT 

Frequency response of an optical modulator, operating at A = 0.83 
pun and utilizing the linear electro-optic ( Pockels ) effect in a Mach- 
Zehnder configuration using an Al g Ga t - t AslAl v Ga,-yAs double het- 
erostructure, is analyzed. We show that in semiconductor modula- 
tors, electroabsorption should be taken into account in optimizing 
the frequency response of the device. C 1993 John Wifev 6 Sons, Inc . 

I. INTRODUCTION 

Light modulators, integrated with optoelectronics, have at- 
tracted intense research activities in recent years. Speed and 
inherent parallelism of all-optical signal processing usually 
justify overall larger size of the optical components, and much 
activity is devoted to establishing fabrication sequences that 
result in reliable integrated devices and waveguides [1, 2], 
Electro-optic modulators based on the linear electro-optic 
effect are expected to have very fast response, limited only 
by the transit times and RCs (R is the resistance and C is the 
capacitance) of the geometry and, in the case of traveling- 
wave modulators, walk off between the light and the micro- 
wave. They are also relatively insensitive to the frequency of 
the light and their RCs and characteristic impedances do not 
change under illumination. On the other hand, the perform- 
ance of electroabsorption modulators is very sensitive to the 
wavelength of the light and their time constants and charac- 
teristic impedances are affected by absorption. However, they 
can be made in much smaller sizes than the electro-optical 
modulators. 

In most cases where one wishes to construct an electro- 
optic modulator from semiconductor materials, some elec- 
troabsorption cannot be avoided. This leads to slower re- 
sponse speeds, because the capacitance of the structure 
increases as electron-hole pairs (EHPs) are optically excited, 
as well as higher than expected attenuation. Particularly, the 
presence of impurity levels in the band gap can considerably 
contribute to the electroabsorption. In this article we analyze 
the performance of a Al^Ga^As/GaAs optical modulator 
based on the above considerations. 


II. THEORETICAL CONSIDERATIONS 

A typical semiconductor-based Mach-Zehnder electro-optic 
modulator is depicted in Figure 1(a) [2]. In this device, the 
integrated optical waveguide separates into two parallel 
branches which form the arms of the interfe rometer. E lec- 
trodes situated in a push-pull configuration produce the elec- 
tric fields which rotate the polarization vector of the light and 
change its phase (the electro-optic effect) in the two arms. 
When the light that is traveling through these two arms is 




S=20pm W=::.5 jjm d= 10 |im ( b) 

Figure 1 (a) Schematic diagram of a Mach-Zehnder modulator, (b) 

The cross section of the Mach-Zehnder modulator fabricated on 
AIGaAs/GaAs 

recombined, due to mode conversion, its intensity is modu- 
lated, 

A cross-sectional schematic through the parallel branches 
[Figure 1(b)] shows the epitaxial structure that constitutes the 
waveguide. It consists of three layers of A^Ga, ,As grown 
on a semi-insulating GaAs substrate. The middle AlGaAs 
layer acts as the guiding channel for the light and the upper 
and lower layers of the AlGaAs act as cladding. The refractive 
indices of different AlGaAs layers are tailored by adjusting 
the molar ratio of the A1 to Ga. The higher concentration 
of aluminum results in lower index of refraction. To minimize 
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the free-carrier absorption, all these layers are undoped. In 
addition, the lower layer also isolates the optical field from 
the GaAs substrate which has a higher refractive index than 
the optical channel. Two-dimensional waveguides are formed 
by etching ridges in the upper cladding layer, which increases 
the effective index of refraction in the region beneath it. The 
epitaxial layer thickness and ridge dimensions, shown in Fig- 
ure 1(b) were chosen to produce single mode guides [3]. Fig- 
ure 2(a) schematically shows ideal electro-optic modulator. 

The induced phase shift in the electro-optic effect is linearlv 
proportional to the electric field (£) and the interaction length 
(Z) [5]: 0 * nlr M EZ (where n {) is the refractive index; r 41 , 
the electro-optic coefficient, for GaAs is 1.1 x 10" 12 m/V). 
Therefore, at photon energies near the band-gap energy of 
the waveguide, the electro-optic_effects are quite small com- 
pared to the electroabsorption which manifests itself through 
the exponential dependence on the absorption coefficient a: 
I(z) — [ { \e~ a - [where I(z) is the transmitted light power at 
position z : /(> is the input power]. Depending on the photon 
energies, a depends on the band structure (band sap), and 
the applied electric field (this dependence has a^ form of 
Ee { ~ CiE) as discussed later). The electroabsorption, then, with 
the exponential dependence on the electric field, can over- 
whelm the electro-optic effect. When electroabsorption oc- 
curs, the absorbed light generates electron-hole pairs (EHP). 
Therefore, the capacitance and the resistance of the structure 
changes under illumination. The change is intensity and fre- 
quency dependent. 

Figure 2(b) shows the band diagram of an n- tvpe semi- 
conductor used in metai-insulator-semiconductor electroab- 
sorptive modulator and Figure 2(c) after the interface. The 
analysis directly applies top-type semiconductors as well. The 
modulating electric field, when it is positive, results in accu- 
mulation at the surface, and when it is negative, it results in 
depletion at the surface. The band bending is exponential in 
accumulation and it has square-root dependence in depletion. 
This yields different probabilities of light absorption, through 
the Franz-Keidysh effect [4], under accumulating and de- 
pleting electric fields. Also, high carrier concentration reduces 
the refractive index, which results in a shift in the confinement 
of the light. The shift is away from the interface under ac- 
cumulation and it is toward the interface under depletion. 

When an insulator is not used between the semiconductor 
and the metal electrode, accumulating electric fields result in 
large current conduction and cannot be used. Thus, in 
Schottky diode structures only depleting electric fields are 
used. This is a serious drawback since a depleting dc bias 
voltage should be applied to prevent forward conduction by 
a microwave. Under reverse bias, dark current flows through 
the diode, and depending on the heterostructure quality and 
the magnitude of the reverse bias, the dark current changes 
as a function of time, deteriorating the modulator perform- 
ance. 

The second problem is generation of photocurrent, which 
increases the power consumption of the device. The photo- 
current is generated due to EHP generation and clearly it is 
light-intensity and frequency dependent. Any impurity levels 
in the band gap will contribute to the current in two ways: 
(i) they increase the dark current, and (ii) they modify the 
carrier lifetime modifying the photocurrent. The most im- 
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Figure 2 (a) A "black box” representation of a light modulator. 

(b) The band structure of a semiconductor under applied electric 
field. Both accumulating and depleting electric field polarities result 
in Franz-Keidysh effect in metal-insuiator-semiconductor structures 

(c) The band-diagram of the modulator. The modulating electric field 
is applied through the high A! concentration layer at the surface 


portant factor determining the response time, however, is the 
capacitance associated with the depletion width in the wave- 
guide and the series resistance. The capacitance is voltage, 
light-intensity, and frequency dependent. As light is ab- 
sorbed, the width of the depletion region becomes smaller, 
increasing the capacitance. 

Using the above considerations, we proceed to determine 
the frequency characteristics of depletion capacitance of a 
Schottky diode when the optical generation rate and trap 
density are nonzero. The electric field dependence of g op can 
be estimated as follows. First, it is easily shown that for unity 
quantum efficiency (i.e., one electron-hole pair per absorbed 
photon) the average number of EHPs generated per unit vol- 
ume per second is given by g ()p = (l/Y)(\I Lhv) fk (/„ - I) 
dx (where V is the volume, L is the interaction length, hv is 
the photon energy, and /’ s are as defined before). Assuming 
a unit area and small a (for below band-gap illumination), it 
is easy to show that g op is proportional to a (g op ~ an/2, n is 
the average number of electron-hole pairs generated per unit 
area per second). According to the Franz-Keidysh effect, the 
electric field dependence of the absorption coefficient is given 
by [5] a = (KE/1E) exp( - C(AE) 1 -/ E) (where K and C are 
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material constants and AE = E K - hi*). Using the numer- 
ical values of the constants in GaAs, and for large fields 
(E/4.74 x 10 4 > 1), we get g op = 1 .24 x 10 -nE. 

Next, using the semiconductor constitutive equations in 
one dimension (.v) in /M ype materials, and using the above 
g op , we arrive at the following well-known linearized fourth- 
order partial differential equation for the potential V inside 
the semiconductor: 

0 4 K/0.t 4 - (1 ItD + iwlD 

+ qiiN D t€ s D) d 2 Vldx 2 - qnafle.D dV/fix = 0, (1) 

where p is the charge density, £$ is the dielectric constant of 
the semiconductor, n is the free carrier density, and 8n is the 
excess carrier density. Moreover, r is the carrier recombi- 
nation lifetime, J x is the current density, q is the elementary 
charge, p is the carrier mobility. D is the diffusion constant, 
E r is the electric field, and N D is the doping concentration 
(all ionized; N D = N5). The direction x is shown in Figure 
1(b). In deriving the above equation, it is also assumed that 
the potential carrier density has a e iuJ time dependence (w is 
the radial frequency of the modulating microwave). Assuming 
that the position dependence of the potential is of the form 
e rx y y’s are given by the roots of the following algebraic equa- 
tion: 

•y 3 — (1 ItD + i<ol D + qfxN D / e s D)y - qi\al2e s D = 0. 

( 2 ) 

In general the potential V is given by 

3 

V(x) = 2 V*v (3) 

i-1 

where V/s are determined from the boundary conditions. 

The following values, which are typical in GaAs at room 
temperature [4], are used in solving for ys: r = 1 x 10 9 
sec, D — 100 cm : /sec. e> = 13 x 8.854 x 10' !4 F/cm, 
N d = 10 i: cm' 3 , AE - 20 meV, q = 1.6c -,19 coulomb, 
k B Tiq = 0.0259 eV (k B is the Boltzmann constant, and T is 
temperature in degrees Kelvin), and D/p = k n T!q. The val- 
ues in AlGaAs, with low concentration of Al. are close to 
the above values, which are more exactly known in GaAs. 
For these values, one of the roots (yj has a positive real part. 
Therefore. V x = 0 for V to stay bounded as.v becomes large. 

The ability of a modulator to properly modulate the light 
is closely related to the ability of an externally applied signal 
to modulate the potential V^.v) inside the semiconductor. In 
Figure 3, y : and y : , as a function of frequency for different 
optical flux density (n) are shown. These plots are equivalent 
to the plots of log( V 2 (x)) and log( V^.v)) as a function of 
log(/) curves at .r = 1 pm. From Figure 3 it can be seen that 
the 3-dB point of the log{ V 2 ) moves from l to 30 GHz as the 
n changes from 10 :1 photons/sec cm 2 10 10 :j photons/sec cm 2 . 
For values of n less than 10 :I photons/sec cm 2 the change in 
the log( V z ) curve is negligible. Of course, this would not be 
the case for different values of r and N D . 

Optical flux density of 10 :: photons/sec cm 2 corresponds 
to a laser with 100-pW power operating at 1.5-eV photon 


energy emitting into an area of Kf ( 1000 x 1000) pm : . With 
n = 10" photons/sec cm~, and the rest of the parameters as 
given above, y versus log(/) plots are generated with r as a 
variable as shown in Figure 4. It is interesting to note that 




Figure 3 (a) y. and y-i as a function of frequency and n (photons/ 
sec cnr). (b) y 2 and y as a function of frequency and carrier lifetime 
r (sec) in GaAs 


shorter r results in wider bandwidths but smaller y> with y 3 
being nearly zero. Smaller y> results in a potential that decays 
much faster inside the semiconductor. This in turn results in 
narrower space charge region where the gradient of the po- 
tential is not negligible. 

The capacitance of the structure is approximately propor- 
tional to the inverse of the space charge width. As r decreases 
or as n increases, the capacitance increases. If the response 
of the modulator is limited by the RC of the structure, larger 
capacitance results in larger RCs and. hence, smaller band- 
widths, 

III. EXPERIMENTAL RESULTS 

The design and fabrication of traveling-wave Mach-Zehnder 
configuration electro-optic modulator, shown in Figure 1(a), 
has been reported elsewhere [1. 2]. Fiber optics were used to 
couple light in and out of the modulator. Optica! coupling 
was accomplished by butt coupling from a single mode fiber 
pigtailed to a laser diode (Ortel LD-620s. peak wavelength 
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of 826.6 mm 1.5 eV). The end facets of the modulator were 
prepared by cleaving; no antireflection coatings were applied. 
The output was monitored by an IR-sensitive camera and/or 
an optical power meter during alignment. 

Measurements with zero applied voltage showed that the 
insertion loss was approximately 24 dB; this includes losses 
from: (1) Fresnel reflection, (2) mode mismatching between 
the fiber and the waveguide, and (3) absorption and scattering 
within the waveguide which includes losses at bends and the 
Y branch. The modulator 3-dB bandwidth was experimentally 



(a) 


(b) 


Figure 4 (a) y ; and y 3 as a function of the optical flux density n 
(photons/sec cm*), (b) y 2 and y 3 as a function of the carrier lifetime 
r (sec) at / = 1 GHz in GaAs 


measured to be 0.5 GHz [2]. The modulator was designed to 
have an 11-GHz bandwidth. 

The observed transmitted optical intensity and the calcu- 
lated cos : <t> intensity, expected for electro-optic modulation, 
as a function of applied voltage is shown in Figure 5. The 
observed intensity versus voltage behavior does not match the 
cos 2 <b behavior, but more closely resembles a continuous 
exponential decrease. Observations on straight sections of 
waveguide (i.e.. without an interferometer configuration) 
showed a similar decrease in intensity with applied voltage. 



Figure 5 Theoretical output intensity versus voltage of a Mach- 
Zehnder modulator and the experimentally measured output intensity 
versus voltage 



Wavefengih ( A ) 

Figure 6 Experimental output intensity versus photon wavelength 
under different applied bias voltages 


Spectroscopy was performed on the waveguides, and Fig- 
ure 6 shows the observed transmission through the waveguide 
as a function of wavelength for applied voltages of 0, 10. and 
20 V. The application of a voltage shifts the absorption edge 
from 814 nm (1.523 eV) at 0 V to 821 nm (1.510 eV) at 20 
V, a shift of 0.013 eV. The transmission drops off at higher 
wavelengths partly due to waveguide cutoff and partly due to 
somewhat lower output of the monochromator and lower de- 
tection efficiency of the detector at larger wavelengths. 

IV. DISCUSSION 

According to Figure 6, the absorption edge at V = 0 is around 
814 nm (1.52 eV). The relationship between the E g and the 
composition index y when 0 < y < 0.35 is [4] E g — 1.424 + 
1.247y (eV). The 1.52-eV absorption edge means that the 
effective y is around 0.08 (8% Al) instead of the expected 
value of 0.1 (10% Al). This discrepancy can also be caused 
by a nonzero internal electric field when V = 0. From the 
position of the laser diode at 826 nm (1.501 eV) it is obvious 
that the applied voltage significantly increases the absorption. 
From Figures 5 and 6 it is clear that the observed intensity 
modulation with applied voltage is due to electroabsorption. 
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Any change related to the electro-optic effect is rendered 
unobservable by the much larger absorption effect. 

To determine the bandwidth of the modulator, we note 
that as the light is absorbed in the Al^Ga^yAs layer, EHPs 
are generated. In the presence of an electric field, EHPs are 
spatially separated. Since this layer has lower A1 concentra- 
tion than the adjacent layers, its band gap is lower. Thus, it 
can contain the EHPs forming a two-dimensional low-resis- 
tivity layer. This sheet of charge provides a relatively good 
conducting layer. Therefore, the effective distance between 
the metallic electrodes is shortened and the capacitance be- 
tween them is increased. Assuming that the top Al r Gai_,As 
layer is semi-insulating, the characteristic impedance of the 
electrodes and the modified structure is calculated, using stan- 
dard microwave programs to be 2.31 Cl. With such a large 
mismatch between the line impedance (50 Q) and the mod- 
ulator. With such a large mismatch between the line imped- 
ance (50 H) and the modulator impedance, most of the mi- 
crowave power is reflected. Only at relatively low frequencies 
is the transmitted microwave able to modulate the potential 
inside the waveguide to any appreciable width. This is dis- 
cussed in the previous section. 

In terms of lumped parameters, the capacitance per unit 
area is 1.15 nf/cm 2 { = eJ2 d, d = 0.5 /xm, and e s = 13 x 
8.854 x 10 14 ), which results in an RC of 1 15 x 10' 1 sec 
(/ *= 0.87 GHz). This is in good agreement with the ex- 
perimentally observed bandwidth of our modulator. 

V. CONCLUSION 

In semiconductors the Franz-Keldysh effect, which results in 
an increase in optical generation rate, can increase the ca- 
pacitance of the structure, lowering the bandwidth of the 
modulator. The effects that are discussed here can also be 
incorporated to increase the figure of merit, for example, by 
including the effect of optical generation rate, and hence the 
excess carrier generation on the refractive index change. This 
also suggests the possibility of using a second light beam to 
modulate light with an applied electric field as an adjustable 
variable that determines the coupling between the two light 
beams through the band bending. It can also be concluded 
that by incorporating an insulator (a semiconductor with a 
large band gap) between the metallic electrode and the chan- 
nel [5, 6], the capacitance of the electroabsorptive modulator 
can be kept low at a reasonable value under the worst con- 
ditions. 
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The dielectric function of a thick layer of Ir^ 52 AIo 48 As lattice matched to InP was measured by 
variable angle spectroscopic ellipsometry in the range 1.9-4. 1 eV. The Ino. 52 AIo. 44 As was 
protected from oxidation using a thin Ino 47 As cap that was mathematically removed for 
the dielectric function estimate. The In 052 Al 04 gAs dielectric function was then verified by 
ellipsometric measurements of other In^ 53 Ga 0 47 As/Iiio 52 Alo. 48 As structures, including 
modulation doped field effect transistors (MODFET), and is shown to provide accurate 
structure layer thicknesses. 


In this letter, we present a measurement of the dielec- 
tric function of In 0 52 A ] 0 4g As and apply it to the ellipso- 
metric characterization of In^Gao 47 As/Ino 52 A 1 0 i 48 As 
heterostructures. The main application of Ino. 52 Alo. 4 gAs is 
as a high band gap semiconductor in In^jGao .47 As/ 
Ino 5 2 Al 048 As heterostructures lattice matched to InP, 
which are used in a variety of microwave and optical ap- 
plications . 1 An important property of Ino. 52 Alo. 4 gAs is the 
fact that, unlike InP, it can be grown in thin film form by 
solid molecular beam epitaxy (MBE), as opposed to phos- 
phorus containing III-V semiconductors which require gas 
sources. Important parameters for any applications that 
can be measured by ellipsometry are the thicknesses of the 
layers, interface quality, and surface contaminations, 
roughness and oxidition. In addition, the dielectric func- 
tion of Ino. 52 Alo. 4 gAs in the visible may be useful in appli- 
cations involving waveguides in this spectral range. 

Ellipsometry, particularly variable-angle spectroscopic 
ellipsometry (VASE) in the visible and near UV, has been 
used to characterize, nondestructive^, a variety of modu- 
lation doped field effect transistors (MODFETs ) 2,3 and 
optoelectronic structures 4 grown on GaAs substrates. Di- 
electric functions of the constituents necessary for the el- 
lipsometric analysis were taken from the literature. How- 
ever, at this time no reliable experimental dielectric 
function of In 0 52 Alo 48 As has been published. Ellipsometry 
has been used twice in the past 5,6 to obtain the dielectric 
function of Ino. 32 Alo. 4 gAs in the visible. In Ref. 5, only 
results for the refractive index were published, while in 
Ref. 6 , the dielectric function of Ino 52 Alo. 48 As was esti- 
mated by scaling the InP values and using an elfective 
medium model with a 3% negative voids fraction. Clearly, 
a direct experimental dielectric function for In^Alo 4S As 
is preferable. 

In most aluminum containing III-V semiconductor 
ternaries, the top layer of the material will oxidize in air 7 in 
a time scale of hours. As ellipsometry is very sensitive to 
the surface conditions, we protected the top layer with 
Ino. 53 Gao. 47 As. We kept the thickness of this cap layer to a 


minimum in order to get reliable results for the 
In 0 52 Alo 4 gAs dielectric function in the near UV, where the 
light penetration depth is very small. As a check of the 
accuracy of our result, we used our experimental 
Ino. 52 Alo. 4 gAs dielectric function to fit two other 
Ino. 53 Gao. 47 As/Ino 52 Alo 4 gAs structures using two parame- 
ter fits only. Finally, we used our result to analyze five 
MBE grown I 1 to. 53 Gao. 47 As/Ino 52 Alo 4 * As complete 

MODFET structures. The MODFETs were grown at two 
laboratories to make sure that no systematic error in the 
analysis was carried over from a systematic error in the 
growth parameters of one group. 

All samples, except three MODFETs, were grown at 
the University of Michigan. The In<). 5 2 Alo 48 As growth tem- 
perature used was the optimum value of 520 *C that was 
found to result in smooth films with the least amount of 
clustering . 8 Some of the MODFET structures were grown 
at 500 "C. The growth rate was in the range 0.6-1. 2 fim per 
hour. The ellipsometric measurements 9 were taken at 3-7 
angles of incidence in the range 300-750 nm (i.e., spectral 
range with reasonable experimental reflectivity) with 10 
nm increments. The calibration sample was measured in 5 
nm increments for better resolution. Marquardt least 
square fits were used to estimate the desired parameters. 

The calibration sample was a U/xm-thick In^Al^As 
on top of a 30 period 3 nm/3 nm Ino . 53 GaQ , 47 As/ 
In 0 5 2 Alo 4 8 As buffer. For ellipsometry purposes, the 
In 0 52 Alo. 4 sAs was treated as the substrate, assuming it to 
be optically thick. The free parameters of the model were 
the oxide and cap layer thicknesses and the values of the 
In 0 5 2 Alo 4 8 As dielectric function at all experimental wave- 
lengths. The calibration functions for the oxide 10 and 
Ino 53 Gao 47 As 11 were taken from the literature. The result- 
ing In aj 2 Al 0 4 gAs dielectric function € is shown in Fig. 1. 
The values of the dielectric function for energies below 1.9 
eV (A > 650 nm) are not as accurate as those above that 
energy due to the light penetrating into the superlattice 
region. Conventionally, a layer is considered optically 
thick if its thickness is greater than 26, where 6 is the light 
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FIG. 1. Dielectric function e of In^AI^As in the range 1. 9-4.1 eV. €, 
and f 2 are the real and imaginary parts of the dielectric function, respec- 
tively. 

penetration depth. 12 We choose to be more conservative 
and use a 3 6 cutoff criterion. The ellipsometric result also 
provided best fit values of 2.2 nm of oxide and 1.3 nm of 
In 0 .s 3 Ga 0> 4 7 As, which are reasonable values for a native 
oxide thickness and a nominal 2 nm In 0 5 3 Gag 47 As cap. We 
obtained a mean square error (MSE) for the tan ¥ and 
cos A fit 3 of 1.0 X 10“ 5 , where ^ and A are the ellipsomet- 
ric experimental results. This exceedingly low value of 
MSE is due to the large number of parameters. 

Next, two samples of thick Ing 52 Alg 4 gAs layers were 
measured, both grown on InP without a superlattice buffer. 
Sample A had a nominal 2-/im-thick Ing^Alg^As while 
sample B had a little over a 1-jmi-thick layer. Both had a 
thicker Ing 33 Gag 4 7 A$ cap layer than the calibration sam- 
ple. The ellipsometric model used for these samples in- 
cluded two parameters only, the thicknesses of the oxide 
and the In^Gag^As cap layer. The Ing 52 Al 048 As dielec- 
tric function used was the result obtained here in the first 
stage of the work. The results for samples A and B are 
summarized in Table I. The values of the MSE for these 
two parameter fits, especially for sample A, are extremely 
good. In both cases, the Ino 53 Ga 047 As layer thickness es- 
timated by ellipsometry is smaller than the nominal value. 
We believe some of the material was oxidized and some 
error may be due to the growth calibration. However, we 
did not encounter this discrepancy in MODFET samples 
made at another laboratory, as will be shown below. 

The MODFET structures shown in Table II were 
made at both the University of Michigan (sample Nos. 1 


TABLE I. Best fits for layer thicknesses in nm for samples made of a 
thick In 052 AI 04 |As layer, considered as substrate. Analysis range 300-620 
nm. 



Oxide 

. I n osjC»ao47As 

MSE 

xio-* 

Sample 

Nominal VASE 

Nominal 

VASE 

A 

2.20 ±0.02 

20 

13.3±0.1 

2.7 

B 

1.3 ±0.1 

40 

32.3 ±0.4 

5.1 


and 2) and by a commercial vendor 13 (sample Nos. 3, 4, 
and 5). All structures had complex buffer layers. For ex- 
ample, the University of Michigan samples had the follow- 
ing layers: starting from the semi-insulating InP substrate, 
a 30 period 3 nm/3 nm Ino 53 Gag 47 As/In 0 ^AI^gAs buffer 
layer and a 400 nm undoped Ing 52 AI 0 48 As as the lower 
part of the conduction channel quantum well. The other 
samples (Nos. 3, 4, and 5) had additional layers below the 
30 period lattice, but they had the same buffer structure 
just below the conduction channel. As the ellipsometric 
analysis was limited to the device active layers, the 400 nm 
In 052 Al 04g As layer was regarded as substrate. Thus, we 
limited our analyses to wavelengths below 540 nm. The 
layer thicknesses, as estimated by RHEED, are given in 
Table II. In all samples, the active layers included an un- 
doped InQ 53 Gao 4 7 As channel, an Ing 32 Alg 4 gAs donor 
layer, and an In 053 Ga 0 4 7 As cap layer. The frig }2 AlQ 4 gAs 
donor layer had a doped (5x 10 18 cm -3 , Si) 15 nm layer 
on top of a 5 nm undoped spacer. The cap layer was also n 
doped, at 3xl0 1R cm~ 3 , Si. In the ellipsometric model, 
doping effects on the dielectric function were neglected. 2 * 3 
The nominal ellipsometric models, including all layer 
thicknesses for the five samples, are given in Table II under 
the heading “Nominal.” A summary of the ellipsometric 
results is given in Table II under the heading “VASE.” The 
errors shown are the 90% confidence limits obtained from 
the least squares fitting. 14 Representative tan ¥ and cos A 
model fits for sample number 5 are given in Figs. 2(a) and 
2(b). In general, the quality of the fits, as given by the 
MSE, is very good, except for sample No. 3. The results for 
the samples made by the commercial vendor are very illu- 
minating. Sample Nos. 3 and 4 were grown in 1991 and 
were nominally equivalent, except for the cap layer thick- 
ness. The VASE results show that in both samples the 
Ing i52 A1 048 As donor layers and the In 0 53 Gag 47 As channels 
are much thicker than the nominal values. However, the 
MSE is much larger for sample No. 3, denoting a poorer 


TABLE II. Best fits for InojjGao^As/Itio j 2 Al a4 |As MODFET layer thicknesses (in nm). Analysis wavelength range 300-405 nm. 4 00 nm 
In 0 s 2 Al 044 As layer used as substrate. 


Sample 

number 

Oxide 

1°0 1.5 3^*0 47 AS 

Cap layer 

Donor layer 

t n 0 530*0 47 AS 

Channel layer 

MSEX10- 4 

Nominal VASE 

Nominal 

VASE 

Nomina! 

VASE 

Nominal 

VASE 

1 

4.6 ±0.1 

10 

4.8 ±0.1 

40 

42.2 ±0.3 

55 

53.9±0.7 

1.6 

2 

2.4 ±0. 1 

15 

8.8 ±0.2 

45 

47.9 ±0.4 

55 

56.1 ± 1.0 

2.4 

3 

2.6 ±0. 1 

5 

2.3 ±0.3 

40 

46.4 ±0.3 

40 

5 1 .2 =±0.6 

4.9 

4 

1.0 ±0.04 

35 

40.6 ±0.4 

40 

48.4±0.5 

40 

54.1 ±1.4 

0.7 

5 

1.4±0.1 

35 

32.9 ±0.5 

40 

43.4 ± 1.1 

40 

4 1.5 ±2.1 

2.4 
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FIG. 2. Experimental and model simulation for (a) tan 'V (b) cos A vs 
wavelength for sample number 5, Table II, for the following three angles 
of incidence; □ 75*, 0 76*. o 77*. 

model. We speculate that the problem in this sample orig- 
inates from the very thin thickness of the cap layer. The 
I n o.53Gao4 7 As top laytfr may be too thin for protecting the 
Ino.52Alo.48As from oxidation, and the sample model shown 
here did not take this into account. Oxidation problems are 
encountered in all thin capped MODFET structures (sam- 
ples No. 1, 2, and 3), where the VASE determined 
Ino 5 3 Gao.4 7 As layer thickness is very small (see Table II). 
The large discrepancy between the nominal and the VASE 
thicknesses of the donor and channel layers in samples 3 
and 4 as compared with that of sample 5 can be explained 
as follows. Sample No. 5 was grown 10 months after the 
other two samples. In the meantime, several improvements 
in the nominal thickness calibrations were implemented, 15 
including repositioning of the RHEED gun, adjusting the 


shutters to reduce transients, and a better correlation be- 
tween RHEED results and ternary alloy thicknesses. In- 
deed, for sample No. 5, the nominal and VASE thicknesses 
are the same to within 8 % for all layers. 

We also analyzed all MODFET samples using the sug- 
gested Ino jjAlo^yAs dielectric function from Ref. 6. The 
results looked unreliable: the MSE were a factor of 5 to 18 
higher for the thin In^Gao^As capped samples; there 
was no consistency between sample Nos. 3 and 4; and most 
thicknesses were far away from the nominal values. 

In summary, we have experimentally determined the 
dielectric function of Ino.52Alo.4gAs in the range 300-650 
nm and have successfully applied it in the determination of 
the active layers’ thicknesses of MODFET devices lattice 
matched to InP grown by MBE. This ellipsometric nonde- 
structive characterization of high performance 
MODFET’s that include Ino ^ jAIq^As layers not only pro- 
vided a confirmation of the nominal layer thickness values, 
but also identified problems in thickness calibrations dur- 
ing growth, as well as cap layer oxidation problems. 
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Abstract 

Metal -organic chemical vapor deposition was used to grow epitaxial layers of AlGaAs, GaAs and InGaAs on 
semi-insulating GaAs substrates. The ternary composition of the thick layers was determined by double-crystal 
X-ray diffraction (DCXRD). Variable angle spectroscopic ellipsometry was used to characterize several types of 
structures including relaxed single-component thick films and strained lattice multilayer structures. The thick film 
characterization included ternary alloy composition as determined by a numerical algorithm and interface quality. 
The results for the alloy composition were equal to the DCXRD results, to within the experimental errors except 
for the top layer of a thick AlGaAs film. The strained layer multistructures were analyzed for all layer thicknesses 
and alloy compositions. For most layers, the thickness was equal to the nominal values, to within the experimental 
errors. However, in all three Iii^Ga^As samples, the indium concentration estimated from the relaxed layer’s 
InGaAs algorithm was around 21%, i.e, much lower than the nominal value. This result indicates a shift in the 
critical points of the dielectric function, owing to strain effects. 


1. Introduction 

The most common technique of growing epitaxial 
III-V semiconductors films on III-V substrates is by 
molecular beam epitaxy (MBE). However, the tech- 
nique of metal organic chemical vapor deposition 
(MOCVD) is more versatile than MBE, as it can be 
used to grow a larger selection of group V materials, 
e.g. phosphorus. In addition, MOCVD can produce 
film growth on several wafers concurrently compared 
with a single wafer in MBE. The materials of interest in 
this study are AlGaAs, InGaAs and GaAs, as they are 
needed to grow the multistructure needed to produce 
modulation-doped field effect transistors (MODFETs). 
This MODFET device was developed by MBE in the 
mid-1980s [1]. MOCVD is being used to grow strained 
layer InGaAs on GaAs [2, 3] but MOCVD growth of 
the combination of a AlGaAs layer adjacent to strained 
InGaAs is not common. The commercial applications 
of this MODFET are widespread [4], especially for 


materials grown on GaAs substrates. The conduction 
channel, made of InGaAs, is under strain. However, if 
the thickness of the channel is below a critical thickness 
h c [5], the strained layer grows pseudomorphically on 
GaAs. These strain effects have essentially no effect on 
the electrical properties of the MODFET, but they can 
be important for optoelectric applications [6], 

The most common technique to calibrate a MOCVD 
reactor for growth rate and alloy composition is 
to measure mechanically or by microscopy the 
layer thickness and to use double-crystal X-ray diffrac- 
tion (DCXRD) for composition. However, DCXRD 
needs a relatively thick layer of material to give a 
reliable value of alloy concentration. A much larger 
thickness is needed to assume a strain-free layer and to 
avoid corrections to the X-ray results, as a result of 
strain. 

In this paper we will show several ways that ellip- 
sometry can be used to help in calibrating the MOCVD 
growth parameters. We will also show instances in 
which ellipsometry was able to pinpoint problems in the 
layer composition or interfaces. In addition, a quantita- 
tive measure of the strain effect on the dielectric func- 
tion of In 03 Ga 07 As will be estimated. 

Thin Solid Films, 
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2. Experimental details 

All samples were grown by MOCVD at Spire Corpo- 
ration, on semi-insulating GaAs(OOl) substrates, using 
low pressure MOCVD in an SPI-MO CVD™450 reac- 
tor. The growth temperature and pressure used for all 
samples were 650 °C and 53 Torr respectively. Source 
regeants used for these layers include trimethylgallium, 
trimethylindium, trimethylaluminum, and 100% arsine 
(AsH 3 ) with a palladium-diffused hydrogen carrier gas. 
Two sets of samples were selected for ellipsometry 
work: (a)thick films; (b) MODFET structures. In set a, 
two aluminum gallium arsenide (Al v Ga,_ v As) films 
were supplied from the same run (M4-2128) with an 
aluminum molar fraction of 18% (averaged over both 
samples) as measured by DCXRD. The Al v Ga l-v As 
film thickness was 8500 A; the growth rate for this layer 
was 2.6 A s' 1 . The Al v Ga,_ v As was protected by a 
thin (50 A) GaAs cap. Two indium gallium arsenide 
(In K Ga t _ r As) compositions were calibrated for Y by 
DCXRD. The first, M4-2093-2 has an indium concen- 
tration of 25.6% and was grown at a rate of 3.5 A s" 1 
for a total thickness of 1.75 pm. M4-21 16-2 is a 1.8 pm 
film of InGaAs with 30% In concentration; the growth 
rate for this composition was slightly higher at 
4.4 As" 1 . Set b contains three samples. All have the 
nominal structure given in Fig. 1, but with slight differ- 
ences. Two of the samples have a 100 A InGaAs chan- 
nel, while one has a 50 A channel. Also, one sample had 
originally a 300 A GaAs cap layer. Some variations in 
the buffer among the samples were also present, but 
their thickness was always 1 pm with at least 9000 A of 
GaAs. As elliposmetry cannot distinguish easily be- 
tween doped and undoped material, we show in Fig. 1 
also the initial model for ellipsometric purposes. 


MODFET STRUCTURE 

Nominal Structure 

Nominal Structure for Ellipsometry 
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Fig. 1. Nominal structure used for ellipsometry compared with the 
actual nominal structure. The ternary concentrations X and Y are 
variables in the ellipsometric analysis. 


The ellipsometric technique was described previously 
[7] and will only briefly be described here. In order to 
increase accuracy, many measurements were made us- 
ing two-zone averaging [8] and/or estimating the angle 
of incidence by using a known sample (GaAs wafer in 
this case) and obtaining the correct angle by least- 
squares analysis. We found these corrections to be very 
small. A model calculation of the ellipsometric parame- 
ters tan *P C , cos A c was least-squares fitted to the exper- 
imental tan *F e , cos using the Marquardt algorithm 
to minimize the mean square error a\ 

= (AT -/>)-'£ [( tan V- e . , - tan f' L , ,) 2 

/ 

-F (cos A e , - cos d e ,) 2 l 

Here N is the number of experimental points and P is 
the number of free parameters in the model. Data at all 
relevant experimental angles of incidence and wave- 
lengths are included. The dielectric functions of 
Al v Ga, _ v As for any value of the aluminum concentra- 
tion X were calculated using the numerical algorithm 
and the critical points given in ref. 10. The functions 
needed for In^Ga^yAs for all Y values were calcu- 
lated using a numerical algorithm similar to that of ref. 
10, but using the following critical points (in electron- 
volts): from ref. 11: E 0 (Y) = 1 .424 - 1 .53 Y + 0.45 Y 2 \ 
E X (Y) was taken as a linear approximation between 
GaAs and the experimental results [12] at Y = 0.3, 
E 2 {Y)—A.%. In this numerical calculation we found 
that very little difference was introduced in the dielec- 
tric functions if the E x critical point was exchanged with 
almost equal energy value of the peak in £,, which is 
much easier to deduce from the experimental result in 
ref. 1 2. 


3. Results and discussion 

All four thick samples were measured as received, 
in the wavelength range 3200-7500 A at five angles 
of incidence each (around 75°). Results for the two 
InGaAs samples are given in Table 1. As these films 
were over 1.75 pm thick, we could not penetrate the 
thickness of the film. Thus we had only two parameters, 
the top oxide film thickness and the indium molar 
concentration Y. In Table 1 we show two results for 
sample M4-2093-2. In the second measurement the 
value of a has improved by a factor to 2.5 because of a 
better angle of incidence definition obtained from mea- 
surement of a GaAs wafer. 

The results for Y are basically equal to the experi- 
mental DCXRD data, to within both techniques’ exper- 
imental error, indicating the reproducibility of our 
results compared with the sample in ref. 12 and the 
validity of the algorithm. The excellent values for o 
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TABLE I. Best fits for Tn ) -Ga ) „ > As thick layers, wavelength range 3200- 7500 A 



M 4-2093-2 


M4-2093-2 


M4-21 16-2 



Nominal 

Best fit 

Nominal 

Best fit 

Nominal 

Best fit 

Oxide thickness 
In concentration (%) 

<7 

25.6 

30.3 + 0.6 
24.1 ± 1.0 
6.4 x IQ" 4 

25.6 

27.9 ±0.3 
24.0 ± 0.5 
2.5 x lO" 4 

30.0 

25.8 ± 0.4 

29.9 ±0.5 
7.8 x 10 4 


denote a smooth film and reliable parameter values. The 
results for the two AlGaAs films are given in Table 2. 
Here f,( A), / 2 (A) and f 3 ( A) are the oxide, GaAs and 
AlGaAs thicknesses respectively and X(%) is the alu- 
minum molar concentration. Comparison of experimen- 
tal and calculated tan *P and cos A for sample M4-2 128-3 
before etching are shown in Fig. 2. An H 3 P0 4 :H 2 02:H 2 0 
(3:1:75) etch was performed on each sample indepen- 
dently. This etch can remove both GaAs and AlGaAs 
material. After the etch, the samples were measured in 
a dry nitrogen environment at an angle of incidence of 
75°, with two-zone averaging. Two measurements were 
made on each sample, one immediately following the 
etch, the other 3 h later. Absolutely no change was 
observed as a function of the measurement delay time. 
Table 2 shows that, with the etch, the top oxide layer 
thickness increased and the AlGaAs thickness decreased 
as expected. We see two problems with the results in Table 
2: the large values of a compared with Table 1 and the 
extremely small values of GaAs thicknesses t 2 before etch. 
In addition, the fits given in Fig. 2 reveal very systematic 
discrepancies between the calculated and the experimen- 
tal tan V' and cos A values in the wavelength range below 
5000 A, i.e. near the surface. 

In order to clarify this last problem, we made fits of 
the data in the range 3200-5500 A. The results are 
shown in Table 3. Here t 4 is the thickness of an interface 
layer of Al 02 Ga 08 As with /,(%) being the volume 
fraction of A1 2 0 3 . The very small amount of A1 2 0 3 
could have been formed either by the aluminum getter- 
ing action during deposition or by post-deposition oxi- 
dation from the air, through the very thin GaAs cap 


layer. The penetration depth of light in Al 0 2 Ga 0K As in 
this range is always below 2000 A. Thus, it turns out 
that the parameter r 3 of Table 2 cannot be estimated 
accurately from results in this restricted wavelength 
range, and therefore the Al v Ga,_ v As was treated as a 
substrate. The most important result from Table 3 is 
that, before the etch, the top layer of the 8500 A thick 
AlGaAs film does not behave like Al 02 Ga OK As, irre- 
spective of the interface model or the quality of fit. It 
seems that a long MOCVD growth of AlGaAs followed 
by GaAs gives a poor quality interface layer, showing 
up in ellipsometry by a larger aluminum concentration 
in the top layer than in the bulk of the film and a small 
amount of A1 2 0 3 . The etch depth of sample M4-2I28-1 
was about 125 A (from Table 2), while that of sample 
M4-2I28-3 was about 225 A. From Table 3 we see that 
the X value is reduced with etching, with higher reduc- 
tion for deeper etch. An approximately 225 A etch was 
required to obtain the same quality AlGaAs at the top 
and bottom interfaces. Then, the ellipsometrically de- 
termined values of X for all wavelength range analyses 
are in excellent agreement with the DCXRD value. This 
result shows that, except for the top approximately 
200 A layer, the AlGaAs film is of good quality, with 
the aluminum concentration constant over all the film 
thickness. We have also performed best fits for the 
wavelength range 5600-7500 A. In these cases we al- 
ways obtained results for X and a similar to those in 
Table 2. The large value of a is probably due to 
inhomogeneities in the very thick AlGaAs film. This 
can be seen from the resonances above 6000 A in Fig. 2, 
where the experimental and the calculated functions 


TABLE 2. Best fits for Al v Ga, _ v As thick films, wavelength range 3200-7500 A 


Sample 

Etch 

a 

i, (A) 

/ 2 (A) 

h (A) 

X <%) 

M4-2I28-1 

Before 

2.32 x 10-’ 

13.8 ±0.7 

0.8 ±2.8 

8515 ± 9 

20.0 ± 0.1 

M4-2128-! 

After 

3.23 x 10-’ 

32.6 ± 1.6 

— 

8392 ±13 

20.7 ±0.1 

M4-2 128-3 

Before 

1.87 x 10-’ 

12.2 ±0.8 

2.5 ±2.5 

8745 ± 9 

19.7 ±0.1 

M4-2128-3 

After 

1.52 x 10-’ 

23.9 ±0.7 

— 

8521 ± 10 

18.1 ±0.2 

Nominal values 




50 

8500 

18 
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Fig. 2. Experimental and model calculated values for (a) tan 'V and (b) cos ,1 tv. wavelength at three angles of incidence for sample M4-2I28-3 
before etching. 


TABLE 3. Best fits for Al t Ga, .,As thick films, wavelength range 3200-5500 A 



M4-2I28-1 Before 1.8 x 10 4 15.8+0.3 

M4-2128-1 Before 7.3x10'’ 18 3 ±0.6 

M4-2128-1 After 6.6x10'’ 3 1 .2 ± 0.3 

M4-2I28-3 Before 2.0 x 10 4 1 3.9 ±0.3 

M4-2I28-3 Before 9.3x10-’ 1 3.6 ± 0.6 

M4-2I28-3 After 8.5x10-’ 22.1 ±0.2 


26.0 ± 1.0 

— 

— - 

2o.o ± 0.4 

13.1 ± 1.4 

78 ±7 

0.8 ±0.3 

30.2 ±0.5 



— 

— 

26.2 ±0.4 

26.3 ± 1.1 

— 

— 

28.1 ±0.5 

14.6 ± 1.6 

75 ±8 

0.7 ±0.4 

29.5 ±0.6 

— 

— 

— 

20.5 ± 0.3 


Nominal values 


50 


18 


1 + 1 + 1 + 1 + 14 - 1 + 


have the same shape, but the magnitudes are a little 
different. 

Results for the three MODFET structures are sum- 
marized in Table 4. Preliminary results were given in 
ref. 9, where no InGaAs algorithm was used. The 
nominal values were obtained from the prior calibra- 
tion of deposition rates. The low aluminum concentra- 
tions in two of the samples are believed to be related to 
a possible problem in these particular runs. For ail 
three samples, we found a correlation parameter a# of 
0.91-0.95 between the GaAs and the AlGaAs thick- 
nesses. For the samples with the lower a fr i.e. M06- 
332-1 and M06-334-1, thicknesses of all layers 
including both InGaAs and AlGaAs are in excellent 
agreement with the nominal values. For sample M06- 
316-1, with = 0.95, the value of the sum of the GaAs 
and AlGaAs thicknesses is within 7% of the nominal 
value. However, in all cases the value of Y, the indium 
concentration in the strained layer, is much smaller 
than the nominal value. For the 100 A strained layer 
samples, we obtain an apparent concentration Y = 23% 
while for the 50 A layer we obtain Y = 18%. This result 
is quantitative, compared with the qualitative results 
obtained in ref. 9 for MBE-grown samples around 
Y = 53%. It must be mentioned that the experimental 
critical thickness h K for In a3 Ga 0 . 7 As is around 90 A. 
Thus the 100 A InGaAs layer may be more relaxed 
than the thinner film, and therefore the strain effect is 
more pronounced for the thinner film. However, as seen 


from Table 4, the errors in Y make the results of all 
three samples overlap, around Y = 21%. 

In order to increase our sensitivity and to decrease 
the 90% confidence limits on the value of Y t we etched 
the thick GaAs layer of sample M06-316-1 and remea- 
sured it at five angles of incidence. An excellent result, 
with <7 = 8 x 1 0 -5 , was obtained using a composite top 
layer, with an oxide layer thickness of 34.7 ± 3.5 A and 
a 72 ± 1 1 A thick layer composed of a mixture of 
AI 0 2 Ga 0 8 As and (11 ±4.5)% of Al 2 O v The remaining 
structure included the AlGaAs layer with f 3 = 
295 ±5.6 A and A" = 19.8 ± 0.4, and the InGaAs with 
t 4 = 94.3 ± 5.6 A and Y = 2 1 .6 ± 0.9. This result for Y 
has the lowest confidence limit. We checked the correla- 
tion of Y with the other parameters, as shown in Table 
5 for sample M06-316-1. No correlation was found, 
except with the InGaAs layer thickness. In all samples 
including the etched sample, the InGaAs thicknesses 
are very near to the nominal values, therefore increas- 
ing our confidence in the values of Y. Thus we believe 
that the value of 7 = 21% for a coherently strained 
In 03 Ga 07 As is reliable. 

This result can be explained by changes in the critical 
points E ( (F), / = 0, 1, 2, with strain e. Recalling the way 
the InGaAs algorithm was structured, we expect E t (e) 
to change in an opposite way to £, (7), i-e. both E 0 (e) 
and £,(£) to increase with £, while E 2 will remain 
constant. The biaxial strain dependence of critical 
points in III-V semiconductors is treated theoretically 


TABLE 4. Best fits for three modulation-doped field effect transistor structure samples, wavelength range 3500-6800 A 


Sample 

a 

i, (A) 

h (A) 

h (A) 

xm 

U (A) 

Y (%) 

M06-316 

2.4 x 10 " 

20.5 ±0.5 

337 ±4 

367 ±6 

19.7 ±0.5 

91 ± 13 

24.1 ±2.7 

M06-316 

Nominal 

— 

300 

450 

20 

100 

30 

M06-332 

5.2 x 10 4 

16.3 ± 0.5 

93 ±3 

437 ±4 

14.3 ±0.4 

94 ±6 

23.5 ±2.1 

M06-332 

Nominal 

— 

100 

450 

20 

100 

30 

M 06- 334 

4.7 x I0~ 4 

16.9 ±0.5 

99 ±3 

488 ±5 

13.7 ±0.4 

65 ± 14 

18.3 ±3.5 

M06-334 

Nominal 

— 

100 

450 

20 

50 

30 


TABLE 5. Best fits for modulation-doped field effect transistor structure sample M06-334 using constant indium concentration values 


Y (%) 

o 

/, (A) 

h (A) 

(A) 

X (%) 

r, (A) 

10 

5.19 x 10- 4 

16.5 ±0.5 

99 ±3 

482 ±5 

13.3 ±0.5 

109 ± 13 

15 

4.74 x I0" 4 

16.6 ±0.5 

99 ±3 

487 ±5 

13.4 ±0.4 

78 ±8 

17 

4.68 x 10 4 

16.8 ±0.5 

99 ±3 

488 ±5 

13.6 ±0.4 

69 ±7 

18 

4.66 x 10 4 

16.9 ±0.5 

99 ±3 

488 ±5 

13.6 ±0.4 

66 ±6 

19 

4.66 x 10 4 

16.9 ±0.5 

99 ±3 

488 ±5 

13.7 ±0.4 

63 ±6 

20 

4.67 x 10' 4 

16.9 ±0.5 

99 ±3 

488 ±5 

1 3.7 ± 0.4 

60 ± 6 

22 

4.70 x 10" 4 

17.0 ±0.5 

99 ±3 

487 ±5 

13.7 ±0.4 

55±5 

25 

4.74 x 10- 4 

17.1 ±0.5 

100 ±3 

485 ±5 

13.8 ±0.4 

49 ±5 

30 

4.88 x 10- 4 

17.3 ±0.5 

100 ±3 

483 ±5 

13.9 ±0.4 

43 ±4 
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in ref. 13. E 0 {e) is split and both E 0 (e) and E,(fi) go up 
with the strain. We found out, by using several values 
for E q in our algorithm, that the calculated value of the 
InGaAs dielectric function in our experimental range is 
only slightly dependent on the value of E 0 (E 0 is below 
U4eV, for Y > 20). Thus, the changes in E () with 
strain can be neglected as a first approximation for 
small values of e„ E, is linear [13] for small £, and so is 
the E, dependence on the indium concentration Y for 
coherently strained layers. Thus the strain dependence 
can be exchanged for composition dependence in our 
algorithm approximation. Therefore the measured 30% 
reduction in the value of Y with strain, together with 
the assumptions of linearity, gives us a simple way to 
measure the correct K c value from the measured ellipso- 
metric result T e , using Y e = 0.7Y c . 


4. Conclusions 

A two-prong ellipsometric study of MOCVD-grown 
layers of AlGaAs and InGaAs was performed, includ- 
ing thick films and strained layer complex structures. 
The study shows that the ternary composition of thick 
non-strained layers can be accurately determined to 
within experimental errors using numerical algorithms. 
In the case of complex structures, thicknesses of all 
layers and the alloy composition of non-strained layers 
can be determined simultaneously, provided that the 
correlations between parameters is no higher than 0.9. 


The composition of strained InGaAs can be estimated 
from the experimental result using a correction factor 
based on a linear approximation of the dependence of 
the critical point E, on composition and strain. 
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Enhancement of Shubnikov-de Haas oscillations by carrier modulation 
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A drastic enhancement of the Shubnikov-de Haas (SdH) pattern is obtained by recording the 
changes in the quantum oscillations of magnetoresistance due to modulation of the carrier 
concentration. The technique enables measurement of the SdH waveform at relatively high 
temperatures and in samples with moderate mobilities. The modulated waveform shows 
selective enhancement of the low-frequency SdH oscillations associated with the upper subband. 

Thus, we were able to recor d very clear oscillations generated by a carrier concentration well 
below 5x 10 10 cm -2 . The theory for this selective enhancement is provided. 


Several fundamental aspects concerning transport 
properties of carriers in a two-dimensional electron gas 
(2DEG) are still unresolved. As the second subband starts 
to be populated, does the electron concentration in the 
ground subband (n ( ) continue to increase? 1 * 2 Is the mobil- 
ity of electrons in the second subband larger than that at 
the ground level? 1 ' 3 Does the carrier mobility drop as soon 
as the upper subband starts filling up 3 or even before? 4 

The primary techniques for determining most trans- 
port parameters are the Shubnikov-de Haas (SdH) oscil- 
latory magnetoresistance, the Hall effect, and conductivity 
measurements. While the total carrier concentration n T is 
obtained from the Hall voltage, the frequency (in 1/2?, 
where B is the magnetic field) of the SdH oscillation ren- 
ders the electron concentration in the subbands of the 
2DEG. The presence of carriers in the second subband is 
manifested by a superposition of the two frequencies cor- 
responding to the concentration of electrons in the ground 
subband n x and in the upper one n 2 . Most frequently n 2 <rt t 
and the determination of n 2 is extremely difficult, in par- 
ticular for lower concentrations {n 2 <5x 10 10 cm' 2 ).* 3 

In this letter, we present a novel technique by which 
we are able to drastically enhance the SdH pattern. The 
technique is based on measuring the change in the magne- 
toresistance due to a carrier generated by a modulated light 
source. In addition to the overall enhancement of the SdH 
waveform, the technique has an important property of se- 
lectively increasing the amplitude of the slower oscillations 
due to carriers in the second subband much more than 
oscillations generated by the ground subband electrons. 

The experimental setup used to generate the modu- 
lated oscillatory magnetoresistance waveform is a conven- 
tional SdH setup with an additional lock-in amplifier and a 
chopped laser. The sample is cooled down in an open-flow 
He cryostat with an optical axis to enable excitation of the 
carriers. The magnetic field is swept up to 1.4 T. The 
change in the longitudinal voltage drop is measured by a 
lock-in amplifier With reference frequency provided by the 
chopper. 

Figure 1 shows the raw data obtained using the mod- 
ulated SdH technique (a), followed by a regular SdH mea- 
surement (b). Both measurements were taken at 4.2 EC on 


an InGaAs/AlGaAs high electron mobility transistor 
(HEMT) structure with a 200 A InQ iGa 09 As well, a 100- 
A-wide Alo^GaoijAs spacer, and a 2x 10 18 cm -3 Si- 
doped barrier of the same composition. The identical fre- 
quency of oscillation of 12.4 T seen in Figs. 1(a) and 1(b) 
indicates a carrier concentration of 6.01 X 10 n cm” 2 in the 
2DEG. The raw modulated waveform is as good as the best 
computed SdH pattern derived after background subtrac- 
tion and normalization. The experimentally modulated re- 
sults are superior to those obtained by taking the numerical 



0.0 0.2 0.4 0.6 0J 1.0 1.2 1.4 


Magnetic field, T 


"’National Research Council-NASA Senior Research Associate, on leave FIG. I. SdH oscillations at 4.2 K on an 1% jGa^^As/Alo ^Ga^ 

from the Dept, of Electrical Engineering and Solid Slate Institute. HEMT structure with n-6.0! \ 10 n cm' 1 and fi,— TT000 cnivV s. (a) 

Tcchnion-lsrael Institute of Technology, Haifa 32000, Israel. Raw modulated data, (h) Regular SdH. 

€> Appl. Phys. Lett. 61 (5), 3 August 1992. Reprinted, with permission, from American Institute of Physics. 
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FIG. X Same ts Fig. I, taken at 9.1 K. 


derivative with respect to the magnetic field, 4 as will be 
discussed below. Figure 2 shows SdH waveforms for the 
same sample measured at 9.1 K. The enhancement of the 
waveform by the carrier modulation technique is obvious. 
The oscillations obtained with this carrier modulation 
method are very clear at this relatively high temperature, 
even though the electron mobility is not too large (p e 
=97 000 craVV s at 4.2 K). 

The technique has an important property of selectively 
increasing the amplitude of the slower oscillations due to 
carriers in the upper subband much more than oscillations 
generated by the ground subband electrons. This selective 
enhancement is demonstrated in Fig. 3 on an AIq jGao 7AS/ 
GaAs HEMT structure grown on semi-insulating GaAs. 
The undoped AJGaAs spacer is 60 A and the dopant con- 
centration is 3X 10 u cm“ 3 . The measurements are taken at 
1.95 K. Figure 3(a) shows the raw data for a regular SdH 
measurement, from which the modulation at the lower fre- 
quency can be extracted only after background subtrac- 
tion. The raw data for modulated SdH is shown in Fig. 
3(b). The amplitude of oscillation at the lower frequency is 
significantly larger than that generated by the ground sub- 
band. The Fourier transform of the modulated SdH shows 
a very clear peak at a frequency corresponding to a second 
subband concentration of 8.4x 10 10 cm“ 2 as compared to 


FIG. 3. SdH oscillations at 195 K on an AlaiGi^As/GaAs HEMT 
structure, (a) Regular SdH. (b) Raw modulated data. Note the clear 
low-frequency oscillations. 


9.7X 10* 1 cm' 2 for the ground one. The Hall concentra- 
tion n r is 1.06X10 12 cm -2 and the Hall mobility is /x f 
= 410 000 cmVv $- 

The selective enhancement can be caused either by an 
uneven distribution of the carriers between the two sub- 
bands, or by the different effect of the excess carriers on the 
amplitudes of oscillation. The ratio between the oscillatory 
resistance A and its zero-field value R 0 for a two sub- 
band system is given by 1 


, * 8 \ Ag, Ag»Ag 2 

— r — =A | \-A 2 ”3 * 

Rq 8q 8q 80 


( 1 ) 


where A x and A 2 are the amplitudes of oscillations at either 
frequency while Bn represents the intermodulation term. 
The zero-field density of states is g$, while is oscillatory 
parts are given by 

f—sir\ (kirn, \ 

Ag,=2« 0 I D T (sX)n P f— J cos • (2 > 

where D r { Y) = K/sii«h( K). X=2irkT/fia) 0 T is the tem- 
perature, io r the cyclotron frequency, r f/ the quantum rc- 
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taxation time, and /’= 1,2 denotes the subband. Practically 
the first Fourier component s = 1 suffices to describe most 
oscillations. 7 

The modulation of the carrier concentration results in 
a change in the oscillatory resistance 5(&R XX ) given by 


(* Rxx \J( dA] 3Bn Ag2 ) — 

\ R 0 J - [\an, + 3/1, g 0 j g 0 


+ 


[>m 


6n , 


f 2 Ag, \ Ag; 

+ |^/Tj + dn 7 g 0 ) g 0 

l A i b **,) f ... 

where the excess concentration 6n T is divided between the 
two subbands as 5/tj and 6n 2 . A straightforward interpre- 
tation for the enhancement of oscillations due to the upper 
subband would be that the electron concentration at the 
ground subband remains almost unchanged once the sec- 
ond subband starts getting populated, and therefore, 
8n 2 >8n l . However, most reported data indicate that the 
electron concentration in the ground subband continues to 
increase, as were our measurements on the samples de- 
scribed above. Moreover, it is reported that 8n { is larger 
than 8n 2 } Therefore, this cannot be the cause of the selec- 
tive enhancement 

Next, the derivative terms in Eq. (3) are investigated. 
The derivative of the density of state is given by 
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This result obtained for carrier modulation should be com- 
pared with the technique of differentiating with respect to 
the magnetic field. 6 Taking the derivative with respect to B 
results in a multiplication by rtf, and since that tech- 
nique greatly enhances the signal due to the ground sub- 
band. It should be emphasized that the carrier modulation 
technique provides experimentally a derivative with re- 
spect to n. Since <o t =qB/m* t both terms in the square 
brackets in Eq. (4) are multiplied \/B. Thus, there is a 
relative increase in the amplitudes at low fields. The first 
term depends on the change in quantum relaxation time 
with subband population. Little data is available on this 
dependence. A theoretical analysis is provided by Isihara 
and Smrcka, 9 from which one can derive that 


TqCC 


(5) 


where € F is the Fermi energy and e r and €„ are the energies 
at which the relaxation time and the conductivity reach 
their maximum. Since the electron concentration depends 


linearly on the Fermi level, Eq. (5) can be transformed to 
a dependence on n. Thus, increases steeply with increas- 
ing population at low values of n since at low energies the 
states become localized and resonant scattering dominates. 
It reaches a maximum at a concentration corresponding to 
and dropping from there on as the electrons can get 
closer to the scattering impurities. 10 Therefore, the contri- 
bution of this derivative term is much larger for the upper 
subband, with its low concentration. The quantum lifetime 
of electrons in the ground subband may have reached its 
peak, i.e., the derivative is close to zero. Thus, the contri- 
bution of this term to the oscillations waveform is much 
larger for the second subband. 

The last part to be analyzed is the derivatives of the 
amplitudes of oscillations and the intermodutation. Expres- 
sions for the amplitudes were derived by Coleridge 1 based 
on the effect of various scattering mechanisms on the re- 
laxation time. It can be shown for the case of n x >n 2 that 
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(6) 


where P u is the intraground subband scattering probability 
and P n is the intersubband scattering probability. It is 
clear from the equation that if A 2 increases, A x decreases 
and vice versa. Following the drop of mobility at the onset 
of population of the second subband, the mobility increases 
with increasing carrier population. Since the dominant 
scattering mechanism is due to the ionized impurities, this 
scattering is reduced by additional carriers screening this 
potential. Thus, it can be deduced that the intrasubband 
scattering decreases with increased carrier concentration. 
Reviewing Eqs. (3) and (6), one can expect a positive 
contribution to the modulated signal generated by carriers 
in the upper subband due to the positive derivative of A 1% 
but a negative contribution to the modulated amplitude by 
excess carriers in the ground subband since the derivative 
of A i is negative. The derivatives of the intermodutation 
are both positive, but as long as the oscillations are not 
very large (as is usually the case when the modulated SdH 
technique is used), this jerm, proportional to the square of 
the oscillations, is negligible. 
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TEMPERATURE INDEPENDENT QUANTUM 
WELL FET WITH DELTA CHANNEL DOPING 

P. G. Young, R. A. Mena, S. A. Alterovitz, 

S. E. Schacham, and E. J. Haugland 

Indexing terms: Field-effect transistors. Transistors , Semicon- 
ductor devices and materials 


A temperature independent device is presented which uses a 
quantum weO structure and delta doping within the channel. 
The device requires a high delta doping concentration within 
the channel to achieve a constant Hall mobility and earner 
concentration across the temperature range 300-1-4 K. Tran- 
sistors were RF tested using on-wafer probing and a constant 

G and F were measured over the temperature range 

300-70 K. 


Fabrication of transistors and Hall bar structures was 
earned out using wet mesa isolation techniques. The tran- 
sistors are a dual, Ti/Au gate design with gale lengths and 
width of 0*8-1 *2 pm and 200pm, respectively. Ohmic contacts 
consist of sequentially evaporated Au/Ge/Au/Ni/Au metal. 

Testing of the devices was carried out to determine material 
transport parameters and transistor performance. The tran- 
sport properties of the structures were measured by Hall and 
Shubnikov-de Haas (SdH) techniques over a tempers ure 
range 300-1 *4 K with a maximum magnetic field of 1*4 T and 
carrier light modulation capability. Transistors were RF 
tested using an HP8510 ANA and an in-house fabricated 
cryostat mounted with Design Techniques probes capable of 
on-wafer probing of the device 5 parameters down to 65 K. 

Results and discussion; Fig. 2 a and b show the measured Hall 
mobility and carrier concentrations for the two devices as a 


Introduction : In most semiconductor devices, the gain changes 
as a function of temperature. For a high electron mobility 
transistor (HEMT), as the temperature is lowered, the gain is 
enhanced due to lower electron scattering while maintaining 
an acceptable carrier concentration. When designing circuits 
and devices for applications requiring operation over a large 
temperature range, this enhancement of the carrier mobility 
has to be taken into account Previous studies on delta doped 
channel structures [1-4] have shown that the delta doping 
creates band bending within the channel allowing for the 
quantisation of the carrier states within the doped region. 

We propose a temperature independent FET with a heavily 
delta doped GaAs quantum well channel using molecular 
beam epitaxy (MBE). In the channel doped structure, tem- 
perature independence is achieved via proper ion sheet den- 
sities within the channel region that degrades the peak 
mobility of the carrier hut creates a constant mobility and 
carrier concentration from 300 to 1*4K. This doping tech- 
nique is contrary to typical HEMT devices where the doping 
occurs outside the channel in order to reduce ion scattering at 
low temperatures. 

Device structures and characterisation : The nominal structures 
grown by MBE for this work are shown in Fig. 1. The con- 
duction channel is a GaAs layer of 150 A confined between 
two Alo.3Ga4.7As layers to form a quantum well. The channel 
layer is delta doped in the centre with silicon to nominal 
concentrations of 1*8 x 10* 2 fan 2 or 6 x 10 I2 /cm 2 . The differ- 
ent dopant concentrations were used to determine the effect of 
the dopant concentration level on the parametric temperature 
dependance. The quantum well structure is used to reduce 
parallel conduction and improve the output conductance of 
the resulting device. 


n* GaAs (350 A) 
undoped AlGaAs (350 A) 
undoped GaAs (75 A) 
undoped GaAs (75 A) 
undoped AlGaAs 
semi-insulating GaAs 


delta doping 
18 x 10 l2 /cm 2 
or 

6 x 10 12 /cm 2 


Rg.1 Device structure cross-section 




Fig. 2 Hall mobility and carrier concentration against temperature 
O mobility 

• carrier concentration 
a Nominal silicon doping 2 x 10 t2 /an 2 
b Nominal silicon doping 6 x 10 I2 /cm 2 
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function of temperature. As can be seen, the Hall mobility for 
the lower doped structure varied by 27% over the temperature 
range with maximum a Hall mobility of 2500cm 3 /Vs. The 
temperature dependence below 125 K is consistent with scat- 
tering from ionised impurities. The carrier concentration 
variation of only 5% is considered to be temperature indepen- 
dent and indicates there is no carrier freeze out 

When the delta doping concentration is increased by a 
factor 3*3 to a nominal level of 6 x 10 12 /cm 2 , the structure 
experiences a fluctuation of the electron mobility of only 7% 
over the same temperature range but at the expense of the 
magnitude of Hall mobility. The Hall mobility was measured 
to be ~2000cm 2 /Vs which is a 20% decrease over that of the 
lower doped sample at room temperature. The carrier concen- 
tration also showed a small 5% variation over the tem- 
perature range. Thus a stable carrier concentration is obtained 
over the temperature range for both dopant levels but only 
the higher doped sample experiences constant HaU mobility. 

Finally, very high 2-D carrier concentrations were observed 
in these structures. For the lower doped sample, a dark Hall 
electron concentration of 2-1 x 10 l 7cm J was measured as 
compared to a value of 2*7 x 10 ia /cro 2 under illumination. 
The large photoconductivity (PC) effect leads to an enhanced 
signal of the SdH oscillation enabling quantitative measure- 
ments and positively confirm the 2-D character of the carriers. 
Quantitatively, we found carrier concentrations under illumi- 
nation for the ground and first excited sub-bands of 
2*1 x 10* Van 1 and 5*3 x 10 u /on J , respectively. This shows 
that for the low doped sample, only two sub-bands are popu- 
lated. We could not resolve all the sub-bands by SdH for the 
higher doped sample, but we expect more sub-bands to be 
populated. This larger number of populated sub-bands is the 
main reason for temperature independent mobility, as the 
electrons in the higher sub-bands are spatially less affected by 
the ionised impurity scattering [5]. 

To demonstrate the temperature independence of the struc- 
ture, FETs made from the highly doped structures were tested 
from 70 to 300 K. The S parameters were measured on- wafer 
to determine F and at 5 GHz as a function of tem- 
perature. and F are shown in Fig 3. An average 

of 7*3 dB was measured over the whole temperature range at 
5 GHz with a variation of only ±0*35 dB. The devices were 



O F w 
t G w 

stable across the frequency range and this constant stability 
appears to be characteristic of the structures with channel 
delta doping. Results for F^ show a similar effect with F„ 
being only affected by system noise and not dependent on 
temperature. 


Conclusion: A device quantum well structure with channel 
doping has been demonstrated with temperature independent 
parameters. Electron HaU mobility and carrier concentration 
were almost constant against temperature independence of the 
RF performance for the transistors. With G ^ and F „ 
remaining unchanged over the whole temperature range, the 
need for temperature compensation is eliminated and cryoge- 
nic design of circuits using the room temperature 5- 
parameters is viable. 

26th May 1992 
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Room-temperature determination of two-dimensional electron gas 
concentration and mobility in heterostructures 
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A technique for determination of room- temperature two-dimensional electron gas (2DEG) 
concentration and mobility in heterostructures is presented Using simultaneous fits of the 
longitudinal and transverse voltages as a function of applied magnetic field, we were able to 
separate the parameters associated with the 2DEG from those of the parallel layer. Comparison 
with the Shubnikov-de Haas data derived from measurements at liquid helium temperatures 
proves that the analysis of the room-temperature data provides an excellent estimate of the 
2DEG concentration. In addition we were able to obtain for the first time the room-temperature 
mobility of the 2DEG, an important parameter to device application. Both results are 
significantly different from those derived from conventional Hall analysis. 


Two of the most important physical parameters gov- 
erning the performance of devices based on two- 
dimensional electron gas (2DEG) are the concentration of 
carriers forming the 2DEG and their mobility. In the past, 
the main effort was placed on increasing the carrier mobil- 
ity by separating the doped layer from the conduction 
channel. This was realized using either homogeneous or 
planar doping in the barrier while the conduction took 
place in the channel, i.e., a modulation doped (MOD) 
structure. Lately the trend is to obtain as high a 2DEG 
concentration as possible. Dopants are introduced into the 
channel, 1,2 drastically reducing the carrier mobility, defeat- 
ing the concept of a MOD structure. There is a tradeoff 
between increased doping level and mobility of the carriers 
in the channel. In all cases, the determination of the 2DEG 
concentration and mobility is an essential part of the char- 
acterization of high speed devices. 

Unfortunately the determination of the concentration 
and mobility of carriers is frequently complicated by the 
presence of a parallel conducting path, in most cases in the 
barrier layer, since it is usually heavily doped. Typical 
room-temperature concentrations in the parallel conduc- 
tion path, AlGaAs for GaAs based devices, and InAlAs for 
InP based devices, are of the order of 10 18 cm* 3 . Even 
though the mobility of these carriers is substantially 
smaller than that of carriers in the 2DEG, in particular in 
MOD structures, their presence may significantly modify 
the Hall coefficient R H and conductivity data. 3 Thus, it 
becomes essential to resort to measurements performed at 
cryogenic temperatures, in order to determine the 2D EG 
concentration, assuming it remains constant as a function 
of temperature. 

In this letter we propose a technique that enables the 
determination of the 2DEG concentration using only 
room-temperature measurements. At the same time the 
room- temperature 2DEG mobility is derived, which is a 
parameter of utmost importance in development of high 
speed devices. The technique is based on recording of both 


^National Research Council -NASA Senior Research Associate, on leave 
from the Dept, of Electrical Engineering and Solid State Institute, 
Techn ion- Israel Institute of Technology, Haifa 32000 Israel. 


the longitudinal and transverse (Hall) voltages across the 
sample versus magnetic field, determining the physical 
magnetoresistance and the change in the Hall voltage as a 
function of field. 

The general expressions for the longitudinal and trans- 
verse resistivities as a function of magnetic field B (in the 
z direction) in the presence of two carriers are given by: 3 

P„=£y/ x =,p 1 j/(£i2+'*u)> (0 

P„=E/J x =R h B= -^ 12 /(D? 2 -M? 2 ), (2) 

with 

1 +o^t?) 

Du — />i + Z) 2 

^12=^1 -M 2 . 

In this letter we assume that the carriers in the conduction 
channel of the heterostructure form a 2DEG and we will 
denote these carriers with index 1 and those in the parallel 
path with index 2. The concentrations of the two carriers 
are n * their effective masses are m ( and the scattering times 
are with 1 = 1,2 for the two carriers. E is the electric field 
while the longitudinal current density is denoted by J x (see 
insert in Fig. 1 ). The Hall concentration and mobility will 
be denoted by /i H and /i H and the Hall scattering factor is 
assumed to be equal to 1. The magnetic field dependence is 
introduced through the cyclotron frequency, (o^—eB/m^ 
where e is the electron charge. At zero magnetic field the 
longitudinal resistivity is equal to the parallel combination 
of the resistivities of the two carriers, while at very large 
fields (<t> f /r,> 1 ) the Hall concentration is equal to the sum 
of the two concentrations. If, for example, the two concen- 
trations are of comparable magnitude, but the mobility of 
the first is much larger than that of the second one, 
the low field Hall concentration equals n v At high fields, 
under the same conditions, the longitudinal resistivity is 
given by: 

Pxx = «2/^|("l+"2) 2 - 
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Magnetic field, T 

FIG. 1. Theoretical fit to experimental data of longitudinal resistivity 
(circles) and Hall coefficient R H (triangles) as a function of magnetic 
field for sample Q12I. Insert shows test configuration. 

In our experimental setup we recorded the longitudinal 
voltage and the Hall voltage continuously versus the mag- 
netic field up to the highest field accessible, 1.4 T. Then, 
using a nonlinear least-squares method we fitted both sets 
of experimental data simultaneously to Eqs. (1) and (2). 
Thus we obtained an estimate of the four variables: n lt n 2 , 
Mi» and /z 2 . Two problems are faced during this process. 
First the Hall voltage increases almost linearly with the 
magnetic field which would result in an improper balance 
of the computation by underweighting the low field and by 
overweighting the high field Hall data. To overcome this 
issue we fitted the Hall coefficient which varies much less 
with magnetic fiekj. A second problem is present in a si- 
multaneous fitting process when one set of numbers is sig- 
nificantly larger than the other. In this case the latter set 
would have a marginal effect on the derived parameters. 
Therefore we introduced a normalization factor. A point in 
mid-field range was chosen for which the resistivity and 
Hall coefficient were calculated. The ratio between these 
figures forms a normalization factor by which the Hall data 
were multiplied, making these data comparable in magni- 
tude to the resistivity data. Upon completion of the com- 
putation the data were renormalized. 

Several aspects regarding the accuracy of the proce- 
dure were examined. The standard deviation of the param- 
eters derived from the fitting process was found to be better 
than 0.5%. Second, the procedure was applied to different 
ranges of magnetic fields to determine the consistency of 
the results. While the scatter between derived parameters 
obtained for different ranges was less than 2% at interme- 
diate temperatures (typically between 50 and 200 K), this 
scatter at higher temperatures may be over 10% in con- 
centration and 5% in mobility. Caution must be taken 
when choosing the range of fit of the Hall data. Since we fit 
the Hall coefficient, which is derived from the measured 
transverse voltage by division with the magnetic field, 
small errors in the low field data may result in a substantial 


TABLE I. Carrier concentrations and mobilities for three samples. Con- 
centrations are in 10* 2 cm -1 and mobilities in cm 2 /V s. 


300 K 4 2 K 


Sample 

n \ 


Mi 

M2 


Mh 

"H 

"»H 

QI2H 

1.09 

3.22 

7570 

1850 

2.80 

5080 

III 

0.96 

QI2I 

1.05 

3.17 

7710 

1920 

2.78 

5150 

1.16 

1.01 

Q11A 

0.557 

4.13 

8080 

1810 

2.92 

4060 

0.624 

0.48 


deviation of the entire fit. Another source of inaccuracy 
may be due to the presence of a large transverse voltage at 
zero field. If the parallel concentration becomes much 
lower than the 2DEG concentration, the accuracy of the 
values derived for n 2 and n 2 may be worse than 10%. This 
is frequently the case with better samples at the lowest 
temperatures (and not exposed to light). Under these con- 
ditions, that the conduction data are dom- 


10x10* 



Temperature, K 


4x101* 



Temperature, K 

FIG. 2. Hall data and parameters derived from fitting procedure as func- 
tion of temperature for sample Q12H: (a) carrier concentrations: (cir- 

cles), n 2 (empty triangles), and n H (full triangles); (b) mobilities: jf, 
(circles), fi 2 (empty triangles), and (full triangles). 
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mated entirely by the 2DEG carriers. The actual accuracy 
of the derived 2DEG concentration was evaluated by com- 
paring it to the results obtained from the SdH and Hall 
experiments performed at liquid helium temperatures. 

Table I summarizes the results of the analysis per- 
formed on three samples. The samples were MBE grown 
MOD structures consistent of AlGaAs barrier with 30% 
A1 and a GaAs well. In samples Q12H and Q12I the bar- 
rier was delta-doped with Si at a concentration of 3.5 X 10* 2 
cm -2 , while in sample Q11A the barrier was homoge- 
neously doped at a concentration of 10 18 cm~ 3 . Sample 
Q12I is identical to sample Q12H except that it was re- 
moved from the GaAs substrate by an epitaxial liftoff pro- 
cess. 4 The 2DEG concentrations are compared with the 
concentrations derived from the Shubnikov-de Haas 
(*s<ih) oscillations and the Hall data (n H ) taken at liquid 
helium temperatures. No second subband population was 
evident from the SdH data on either sample. 

For all the samples examined at room temperature, the 
derived 2DEG concentration rtj falls between the measured 
low-temperature SdH and Hall concentration. The differ- 
ence between these and the room-temperature Hall con- 
centration is very significant. Altogether the high- 
temperature Hall data are much closer to that of the 
parasitic parallel layer. It should be pointed out that this 
large parallel concentration n 2 is typical for these types of 
structures. The concentration n 2 drops drastically with 
temperature and reaches 2X 10 M cm' 2 below 50 K. The 
relatively larger error in room temperature n x for sample 
Q12H is due to a large zero field transverse voltage (0.1 
mV with a current of 2 /iA). At the lowest field this volt- 
age is more than half of the recorded signal. 

Figure 1 shows the experimental points and the results 


of the fitting process to the data measured at room tem- 
perature on sample Q12I. The fit of p** is excellent, while 
that of /2 h exhibits some fluctuations at the lowermost 
magnetic fields, due to the division by B , as explained be- 
fore. Even for these points, the error is below 2%. 

As an extension of the method we used it to estimate 
the four parameters n 2 , and /i 2 as a function of 
temperature. Figure 2(a) shows the carrier concentrations 
and Fig. 2(b) the mobilities as function of temperature for 
sample Q12H. The low-temperature 2DEG parameters al- 
most coincide with the Hall data because of the low par- 
allel concentration in this temperature range. It is interest- 
ing to note that the mobility of the parallel layer p 2 
increases as n 2 decreases with descending temperatures. 
However, \i 2 remains fairly constant at lower temperatures 
as might be expected in degenerate semiconductors. 

In conclusion we have shown that the simultaneous fit 
of the magnetic field dependent longitudinal and transverse 
resistivities can serve as a useful tool in characterization of 
quantum well structures both enabling room-temperature 
determination of 2DEG concentration and mobility as well 
as exploration of the effect of parallel conductance. This 
may be particularly helpful in better understanding of per- 
sistent photoconductivity in these structures. 5 
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